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ANALYSIS OF LIQUID-METAL-JET IMPINGEMENT COOLING 
IN A CORNER REGION AND FOR A ROW OF JETS 
by Robert Siegel 
Lewis Research Center 

SUMMARY 

Jet impingement is a possible means for providing effective localized heat transfer. 
Two configurations are analyzed here for jet impingement against uniformly heated sur- 
faces. The first is for flow into the interior or over the exterior of a corner formed at 
the apex of two flat plates. The second is for an infinite row of parallel slot jets imping- 
ing normally against ari)late. Liquid metals are very effective heat -transfer fluids and 
are specifically studied here as a possible means for obtaining high heat -transfer coeffi- 
cients. An analytical aid in the study of liquid -metal heat transfer is the fact that the 
Prandtl number is small, and hence viscous diffusion effects are small compared with 
the molecular diffusion of heat. To a good approximation the flow can then be assumed 
inviscid with regard to the solution of the energy equation. For inviscid flow the jet re- 
gion can be analyzed by conformal mappii^. The flow region is transformed into a re- 
gion of uniform flow between parallel plates. The energy equation is also transformed 
into the coordinates of the parallel plate region and takes the form of an equation for 
which the exact solution is known. The solution is then transformed back into the phys- 
ical plane to obtain the local heat -transfer characteristics along the impingement sur- 
faces. The heat -transfer behavior is given for incident jet Peclet numbers from 5 to 50. 
The results show how the heat transfer improves as the included corner angle or the jet 
spacing is increased. 


INTRODUCTION 

Jet impingement is a possible means for obtaining effective localized coolii^. The 
specific type of fluid considered here is a liquid metal, which is a very effective heat- 
transfer medium that can be used at high temperatures without requiring high pressures. 
In this report some effects of impingement plate and jet geometry are investigated by 
analyzing the heat transfer for impinging slot jets in two configurations. One of these is 
the somewhat confined flow geometry of a single slot jet impinging in a corner region 



formed by two flat plates joined at an angle to each other. The heat -transfer behavior is 
obtained as a function of corner angle. Since the same solution applies for both flow into 
the interior of the corner and over the exterior of the corner, results were computed for 
a wide range of corner angles (6 = 60° to 360°). The second geometry considered is an 
infinite row of uniformly spaced parallel slot jets impinging against a flat plate. It is 
desired to obtain the heat-transfer dependence on the spacii^ between the jets. This 
geometry was studied for air in reference 1, 

The present work is an extension of reference 2 where the impingement heat trans- 
fer was obtained for a single jet flowing normally against a flat plate. The solution in 
reference 2 is a limitii^ case for both geometries treated here; for the corner flow it 
corresponds to an included angle 5 of 180°, while for the row of parallel jets it corre- 
sponds to infinite spacing between the jets. The theory is based on the use of conformal 
mapping, which is valid for inviscid flow; hence, the theory utilizes a potential flow ap- 
proximation with regard to the fluid behavior. Inviscid flow is a good approximation for 
liquid metals for the jet Peclet number range treated here (Pe < 50). These fluids have 
very small Prandtl numbers (Pr 0. 005 to 0. 02). Because the Prandtl number is the 
ratio of viscous to thermal diffusivity, for a situation where both the viscous and ther- 
mal boundary layers are developing simultaneously, the viscous layer will be very thin 
compared with the thermal layer thickness. The heat transfer can thus be computed 
quite well by neglecting the viscous boundary layer and using an inviscid approximation 
with regard to the flow. 

The free streamlines bounding the jet, and the fluid velocity alor^ the heated bound- 
ary (this velocity is needed to evaluate the solution of the energy equation) are obtained 
by mapping the jet region into a potential plane. In this plane the flow geometry be- 
comes a strip of constant width. The mappir^ is achieved by use of an intermediate 
hodograph (complex velocity) plane. 

The strip region into which the jet maps in the potential plane is a convenient region 
in which to solve the energy equation. When the energy equation is transformed into the 
coordinates of the potential plane, it has the same form as the energy equation for a 
\miform velocity in a parallel plate channel. The bovindary condition of uniform heating 
at the impingement plate transforms into a nonuniform heating condition along one of the 
channel boundaries. The analytical solution is available for this situation, so the solu- 
tion is obtained in the potential plane and then mapped back into the physical geometry. 
This yields the temperature distribution aloi^ the uniformly heated wall. The heat 
transfer performance at the stagnation point is also expressed as a heat -transfer coeffi- 
cient, and its behavior is compared with that for single -jet flat-plate impingement. 

This clearly shows the effect of corner angle and spacing between jets in an infinite row. 
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ANALYSIS 


Description of Physical Situation 

The flow configurations, which are two-dimensional, are shown in figures 1 and 2. 
(See appendix A for symbol definitions. ) hi figure 1 a slot jet having a half -width b im- 
pinges against a wedge-shaped plate either into the corner or over the outside of the cor- 
ner, the special intermediate case being impingement against a flat plate. The plates 
forming the corner have a uniform heat flux supplied along them. There is a stag- 
nation point where the centerline of the jet meets the apex of the corner, and it is the 
heat-transfer behavior in this region that is of most interest. The impingement plate is 
assumed to extend to infinity in the direction of the flow along it, but any plate length 
that is longer than several jet widths would provide the same flow pattern in the stagna- 
tion region. The nozzle forming the jet is far enoi:gh from the plate so that the impinge- 
ment does not alter the uniform flow leaving the nozzle. For impingement of a single jet 
against a flat plate, this is true for nozzles more than about three jet widths away from 
the surface. After turning, the flow moves out along each half of the plate, and, for the 
condition of inviscid-irrotational flow considered here, the flow becomes parallel to the 
plate with an asymptotic thickness b. 

In figure 2 the geometry is shown for an infinite row of parallel slot jets striking a 
uniformly heated wall. The half width of the incoming jets is b, and there is a imiform 
half spacing S between them. 


Flow Characteristics of Jet 

The fluid being considered here is a liquid metal, and hence it has a very low 
Prandtl number (Pr « 0. 005 to 0. 02). The Prandtl number is the ratio of the fluid vis- 
cous diffusivity to the molecular diffusivity of heat. In a boundary -layer type of flow 
over a heated surface where both the viscous and thermal boundary layers are developing 
simultaneously, the viscous layer in a liquid metal will be very thin compared with the 
thermal layer thickness and can be neglected. A flow that is initially irrotational, as is 
the case for the approaching jet, will develop little vorticity and can be assumed to re- 
main irrotational. Consequently, the jet configuration can be obtained from free stream- 
line theory by use of conformal mapping. The details of the mappings are given in 
appendixes B and C, and only the final results are summarized here. 

Corner flow . - Because of symmetry, only half of the jet need be considered. It is 
shown with some additional notation in figure 3(a) and in dimensionless form in fig- 
ure 3(b). In the mapping the flow region is transformed from the dimensionless physical 
plane into a potential plane in which the free streamline 12 and the streamline 345 along 


3 



the jet centerline and the heated plate become parallel lines (fig. 4). Thus in the poten- 
tial plane the flow region is a parallel plate channel where heating is being applied to a 
portion of one side. As will be given later, the energy equation is solved in this channel 
flow region, and to transform the calculated plate temperatures from the potential plane 
to the physical plane the relation is needed between $ in figure 4 and L in figure 3(b). 
This relation is obtained from the mapping in appendix B as 


-4 

*'0 




^^(1 - 


0 ^ ^ ^ 1 


(la) 


$ = — tanh“^(4^) 0 4 — 1 


(lb) 


Tf 


To obtain an indication of the flow pattern, the free streamlines were calculated. 
The coordinates are given by 



cos 

0(1 + /S) - /3 I 
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sin 2e 


de 0 < 0 < 1 
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(2a) 



sin 0(1 + 0) 




sin 20 


d0 0^0:^- (2b) 
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The Xq and are given by (these are the X and Y at 0 = »/4) 


X^=icos + 


F(^) 

(3a) 

F(I3) 

(3b) 
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where 


F(/3) = 


Iq 


Row of jets. - From symmetry it is only necessary to consider the region between 
the centerlines of a jet and an adjacent backflow region (fig. 5). In the potential plane 
this region maps into a strip in the same fashion as the corner flow discussed pre - 
viously. The heated wall 71 occupies a portion of the zero streamline. The mapping 
between figures 5 and 6 is given in appendix C and is in terms of elliptic functions. The 
relation between X and $ along the heated wall is 


where 


= - [tan~^ I ■) + 

7 T L V cn £ / J 2b 2b 


$=iln (1 + - lln /I ^ 

7T L(l-k)sn4J ff \1 - k/ 7T \l-k/ 


-K < 4 < K 


To obtain k and k*, the K is first obtained from the jet spacing using the relation 


l/A . i\ 1 

2 \2b / 2 
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The k is then found from K by using the complete elliptic integral 
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and k* is related to k by 


k’= Vi -k^ 
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The coordinates of the free streamlines are given by 


and 


where 



(8a) 


(8b) 


The Governing Energy Equation and Boundary Conditions 

The energy equation in the two-dimensional flow region is given as a balance of 
convection and conduction by 

pCpU • Vt = K V^t (9) 

For inviscid-irrotational flow the fluid velocity can be expressed as the gradient of a 
potential (p: 

u = V<p (10) 

where <p satisfies Laplace’s equation 

V^<p=0 (11) 

The velocity components are related to the stream function by the Cauchy -Riemann 
equations: 

u = (I2a) 

9x 3y 
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and 


y = (12b) 

9y 3x 

By use of equation (10) the energy equation then becomes 

pCp Vcp- Vt = K V^t (13) 

The thermal boundary conditions are as follows (the boundaries without parentheses 
correspond to fig. 3(b) while those in parentheses correspond to fig, 5). Aloi^ the cross 
section of the incoming undisturbed jet, the fluid is at a uniform temperature t^; 

t = *oo x,y on 23 (56) (14) 

It is assumed that along the free streamline there is little heat loss in comparison with 
the convective heat flow within the jet, so this boundary is assumed insulated: 

iig • Vt = 0 x,y on 12 (345) (15) 

Boundaries formed by symmetry lines have no heat flow across them and hence also act 
as insulated boimdaries: 


^ = 0 x,y on 34 (12 and ^) (16) 

9x 


Along the solid boundary a uniform heat flux q^ is imposed so that 


= - — x,y on 45 (71) 


an 


w 


(17) 


Before obtaining the solution, the energy equation and boundary conditions are 
placed in dimensionless form, Usii^ the dimensionless variables defined in appendix A 
the energy equation becomes 


V# • VT = V'^T (18) 

2 


with the boimdary conditions 
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X, y on 12 (345) 


(19) 


A*. 



VT 


= 0 


^ = 0 x,y on 34 (12 and ^) 

ax 


( 20 ) 


aT 

aN 


= -1 


X, y on 4 5 (7 1) 


w 


( 21 ) 


The strip regions in figures 4 and 6 have a simple geometry and hence are conven- 
ient regions in which to solve the energy equation. Since in the potential plane stream- 
lines with equal numerical increments are spaced equal distances apart, the strip rep- 
resents a channel havir^ a uniform velocity distribution across its width. The energy 
equation (eq. (18)) has the same form as that of equation (16) of reference 3, where it is 
shown how the equation can be transformed from the dimensionless X, Y plane to the 
coordinates. Using equation (26) of reference 3 yields the energy equation as 


Z? IT 
2 0$ 


0^T 0^T 

04-^ 


( 22 ) 


This is the same as the energy equation that applies for forced convection in a uniform 
velocity flow in a channel between parallel plates. The flow is in the 4> direction, and 
the channel width extends across the 4^ direction. If the solution is obtained in the po- 
tential plane, it can then be mapped into the physical plane since the $,4^ plane is re- 
lated to the X, Y plane by conformal mappir^. To solve equation (22), the boundary 
conditions in the physical plane must be transformed into the 4>,'If coordinate system. 

The free streamline along which equation (19) applies is a line of constant and 
the normal direction to <!f would thus be along a constant 4> line. Hence, the boundary 
condition in equation (19) becomes 

= 0 4.,^ on 12 (345) (23) 

04 ^ 

The symmetry condition in equation (20) is preserved under the transformation to the 
potential plane. This boundary condition then becomes 

^=0 4,^ on 34 (12 and ^) (24) 

04 ^ 
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To transform equation (21), the fact that the normal coordinate is a function of $ 
and ^ is used to yield 


9T _ 9T + 

3N^ 3N^ a# 3N^ 

Since the derivative of $ in a direction is the velocity in that direction, the 8#/ 8N^ 
at the wall would be zero as there is no velocity normal to the wall at the wall. Then 
along 45 using equation (2 1) 


8^ 8N^ 


w 


K L and form a coordinate system tangential and normal to the wall, the Cauchy- 
Riemann equations yield 


8 ^ _ ^ _ y 

^~8L L 

where Uj^ is the tangential velocity along the wall. The boundary condition then 
becomes 


H = -J_ on45(^) (25) 

^ Ul 


From the solutions in appendixes B and C the Uj^ for the two geometries are given by 



on 45 in figure 4 


$ on 7 1 in figure 6 


(26) 


(27) 


The jet approaching the heated plate is at a imiform temperature, so in the transformed 
plane 

T = T^ 4>,'F on^(^) (28) 
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Solution for Wall Temperature 


Considering the situation in the potential plane, the term on the left side of equa- 
tion (22) represents convection along the axial (that is, #) direction in figures 4 and 6, 
and the terms on the right side are axial and transverse heat conduction. In figures 4 
and 6 the geometry is a moving slab between parallel plates. The slab is at uniform 
temperature for large negative $, and the parallel planes are insulated except for a 
nonuniform heat addition aloi^ the positive 4>-axis in figure 4 and along a portion of the 
$-axis in figure 6. The solution for this situation including axial heat conduction was 
given in reference 2. The specific quantity of interest is the temperature along the wall 
45 in figure 3 or 71 in figure 5, and this is given by 



where for the case in figure 3 


Tjj = 0 and 772“°° 


and for the case in figure 5 
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Except at the value tjg = °°» zero at the r}^ and 7^2 as these correspond 

to the jet stagnation points or central point on the plate of the recirculation region. Al- 
though the Integrands become infinite, the singularities are integrable. For example, 
for small ?j the in equation (26) becomes 



yi+i8)/2 


= ^(l+ 3)/2 


1} « 1 


Then with tJj = 0 the first integral in equation (29) can be written by using a small € 

/»$ 


A'*’ 

_ / d-n , / dr} _ 2 

^r(v) I „(l+^)/2 / Ul(7?) 1-^ 

. 0 ^ ‘^77=0’’ «^e ^ 


2 ^( l -^)/2 ^ 


f ^ 

1 "Lfl* 


This form shows that there is no starting difficulty in the integration and that the correct 
solution can be obtained by starting a small e away from the limit where Uj^ goes to 
zero. This type of behavior is true at all the limits where ~0. For example for the 
infinite row of jets the first integral becomes 


where 





drj 

DlW 



Wall temperatures for large values of L. - For the geometry in figure 3 a solu- 
tion of simple form can be obtained for large values of L. This is done in appendix D. 

In figure 20, which is used in that analysis, it is noted that for large L (where L = £/b) 
the inviscid jet flow becomes a uniform flow with thickness b. The free streamline 
boundary is assumed insulated, and the wall has a uniform heat input. Thus for large L 
the system can be analyzed as a channel flow with uniform heat input. A heat balance 
starting from L = 0 determines the mean temperature at any L. The resulting wall 
temperature for large L is . 
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( 30 ) 


T(L)-T =-?_L+i 
” Pe 3 

It is noted that this is independent of the angle /3 of the plate. 

Local Nusselt n umber. - The local heat-transfer behavior can also be expressed as 
a local heat-transfer coefficient or Nusselt number. The local heat-transfer coefficient 
is 


h = 


t - 1 

W oo 


(31) 


Then the local Nusselt number based on the jet width 2b is 


Nu(L) = 


h2b 


q 2b 


VL)-T, 


(32) 


RESULTS AND DISCUSSION 
Jet Flow Into or Over a Corner 

The free streamlines of the jet were calculated from equations (2) and (3) for various 
angles of the impingement plate. (The plate angle y or corner angle 6 are defined in 
fig, 1. ) Results are given in figure 7 for positive and negative plate aisles. The case 
y = 0 corresponds to impii^ement from the normal direction against a flat plate. For 
y = -tt/ 2 the heated wall is parallel to the jet and is alor^ the jet centerline. In this in- 
stance there is no impingement as the flow direction is parallel to the wall. Figure 7(a) 
shows that the penetration of the flow into the corner region between the plates becomes 
poor as the included angle 6 between the plates is decreased. This will yield decreased 
heat -transfer performance for small values of 6 (large positive values of y). 

Figure 8 shows the waU temperature variation computed from equation (29) for im- 
pingement against plates at various angles. Since the plate is uniformly heated, a region 
of poor heat transfer in the vicinity of the stagnation point 3 delds a high wall temperature 
in that region. As the angle y is increased, there is a rapid increase in temperature at 
the stagnation point £ = 0. This behavior would be expected from the free streamlines 
in figime 7 (a). For y = 60° the large distance between the free streamline and the 
origin indicates a very low flow velocity in the corner region. The parts (a) to (d) of fig- 
xme 8 correspond to the jet Peclet munbers, 5, 10, 20, and 50. An increase in Peclet 
number can be regarded as an increase in jet velocity. An increased velocity produces 
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improved heat transfer; hence, the plate temperature decreases as the Peclet number 
increases. As shown by the final result in appendix D the dimensionless temperature 
for large L depends only on L and the Peclet number. Hence, all the curves for each 
part of figure 8 go toward the same limitii^ line as the L is increased and this line ap- 
proaches the line predicted in appendix C. 

The results shown in figure 9 are a continuation of those in figure 8. In this in- 
stance the plate armies y are negative (6 > 180°) so that the flow is over the outside of 
the heated corner. The heat transfer improves continuously as the plate angle becomes 
more negative, with the best heat transfer being when the flow is parallel to the heated 
plates, that is, when y = -90°. 

The heat-transfer characteristics at the stagnation point are examined more com- 
pletely in figures 10 and 11, Figure 10 gives the dimensionless wall temperatures at the 
stagnation point as a function of plate angle, each cixrve being for a different Peclet num- 
ber. For a fixed Pe the wall temperature increases continuously with y as a result of 
the decreasing local velocities in the vicinity of the stagnation region. 

As given by equation (32), the Nusselt number is proportional to the reciprocal of 
the dimensionless wall temperature. Figure 11 gives the ratio of the heat -transfer co- 
efficient at the stagnation point for flow into or over a corner as compared with that 
for impingement against a flat plate. This shows quite well the significant decrease in 
heat-transfer coefficient as the geometry varies from flow parallel to a flat plate 
(y = -90°) to a jet impinging into a vee groove with total internal angle 6 of 60° 

(y = 60°). For the case of flow parallel to the plate (y = -90°), the usual boimdary -layer 
solution for the energy equation gives an infinite heat-transfer coefficient at the leading 
edge. This is a consequence of neglectii^ axial heat conduction, which means that the 
fluid contacting the leading edge is at the free-stream temperature and the thermal 
boundary layer at the leading edge is of zero thickness. The present analysis however, 
includes axial conduction in the liquid metal so that some heat is conducted upstream in- 
to the fluid before its arrival at the leading edge. This yields a finite heat-transfer co- 
efficient at the leading edge, L = 0. 


Infinite Row of Impinging Jets 

Figure 12 shows the jet free streamlines for several spacings between adjacent slot 
jets in an infinite row. As shown by figure 2 the region in figure 12 is between the two 
centerlines for an incident jet and an adjacent recirculation region. These centerlines 
are at X = ±S/2b and are shown as an example in figure 12 for S/b = 5. 2; to simplify 
the figure the centerlines have been omitted for the other S/b values. As shown by the 
analysis in appendix C , the free streamline at X = 0 moves upward to an infinite Y 
value as S/b approaches 4 so that within the constraint of potential flow for S/b -*-4 
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the region between the jet and adjacent recirculation region becomes completely filled 
with fluid. At large spacings the recirculation region will not affect the incoming jet 
flow, (This condition has been reached to a good approximation when S/b = 8, ) 

The wall temperatiire distributions for impingement on a uniformly heated plate are 
shown in figure 13 for various spacings and jet Peclet numbers. The best heat transfer 
is at the stagnation point of the incident jet (f = 0 in fig, 2 or L = 0 in fig, 5), and poor 
heat transfer is obtained at the center of the recirculation region. The wall temperature 
at the center of the recirculation region goes through a minimum as the spacing is in- 
creased, This results from the improved flow distribution as indicated by figure 12; for 
larger spacings this is coimteracted by the greater extent of the heated region, and the 
temperature rises with spacing. The results for a single jet, which correspond to an in- 
finite spacing, are from reference 2, For a spacing S/b of 6, 4 the results in the vi- 
cinity of the incident stagnation region are within a few percent of the values for infinite 
spacing. 

Figure 14 shows the Nusselt number at the stagnation point of the incident jet (L = 0 
in fig, 5), As shown by the free streamlines in figure 12 the jets significantly interfere 
with each other for S/b less than about 5, 2, causing a thickening of the flow region and 
hence lower velocities in the vicinity of the heated plate. Figure 14 shows how the heat 
transfer becomes poorer as this interference increases with decreased jet spacing. As 
S/b approaches four, the Nusselt number decreases toward zero as the region between 
the jet and adjacent recirculation region becomes filled with fluid. 


CONCLUSIONS 

The flow in an impinging two-dimensional free jet (a single jet or one jet of an in- 
finite row) is bounded by free streamlines and solid boundaries or symmetry lines which 
are also streamlines. As a result of being bounded by streamlines, the jet region maps 
into a parallel flow region in a system using streamline coordinates. When the energy 
equation is transformed into the streamline coordinate system, it becomes the same as 
the equation for convection to a uniform flow in a parallel plate channel. This equation 
can be solved by known analytical techniques, and the solution is then mapped back into 
the physical plane. Thus the analysis combines the conformal mapping technique for 
free jets with a transformation of the energy equation to reduce the problem to one hav- 
ing a standard type of solution. 

Heat transfer results are obtained for two geometries. The first is for a jet im- 
pinging either into a corner region formed at the apex of two uniformly heated plates at 
an angle to each other or flowing over the exterior of the corner region. The second is 
for an infinite row of parallel jets impinging against a uniformly heated plate. 
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The heat transfer for flow into a corner becomes quite poor as the Included angle 
between the heated plates is decreased because a large region of low velocity is formed. 
Results are given to illustrate the continuous decrease in heat -transfer performance as 
the geometry is changed from flow over the outside of a corner (included angle > 180°) 
to flow against a flat plate (included angle = 180°) and then into an internal corner (in- 
cluded angle < 180°). 

In a similar fashion the heat-transfer performance decreases for a row of impinging 
jets as the spacing between jets decreases. When the ratio of spacing between jets to 
jet width increases beyond about 6, the heat-transfer behavior near the impingement 
stagnation point becomes close to that for a single jet. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 1, 1975, 

505-04. 



APPENDIX A 


SYMBOLS 


b 

s 

cn,dn, sn 

f 

h 

A A 

ij 

K 

K' 

k 

k’ 

L 

£ 

m 

N 

Nu 

n 

A 

n 

Pe 

Pr 

Qw 

Re 

S 

T 

t 

t 


half width of incident jet far from plate 
specific heat of fluid 
constant evaluated in eq. (D6) 
constants of integration 
elliptic functions 
function of coordinate n^ 
local heat-transfer coefficient along plate 
unit vector in x and y directions 
complete elliptic integral of first kind 
complete elliptic integral K* = K(k’) 
modulus of elliptic function 
complementary modulus, k’ = f^l - k^ 
dimensionless coordinate, £/b 

coordinate aloi^ heated plate measured from impir^ement stagnation point 
integer 

dimensionless normal coordinate, n/b 

Nusselt number, h2b//c 

normal coordinate 

unit normal vector 

jet Peclet number, |v^|2b/o! 

fluid Prandtl number, C p/k 

heat flux specified at wall 

jet Reynolds number, | v^ 1 2bp/ p 

half spacing between jets 

dimensionless temperature, t/c/bq^; variable in auxiliary T -plane 
temperature; variable in auxiliary t-plane, t = ^ - iri 
mixed mean fluid temperature 
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U, V dimensionless velocities, u/|v^|, v/|v^| 
u velocity in x direction; variable in auxiliary u-plane 

V velocity in y direction 

u local fluid velocity vector 

W complex potential W = # + i^ 

X,Y dimensionless coordinates, x/b and y/b 

x,y rectangular coordinates 

Z complex variable, Z = X + iY 
a thermal diffuslvity of fluid, /c/pC 

Jr 

/3 plate angle ^ = r/{n/2) 

y plate angle measured from direction normal to jet, see fig. 1 

6 total included corner angle (defined in fig. 1) 

€ a small number 

^ complex velocity, ^ = U - iV 

j] dummy variable of integration; ordinate in auxiliary t-plane 

e angle in auxiliary t-plane (defined in fig. 16) 

K thermal conductivity of fluid 

fjL fluid viscosity 

V fluid kinematic viscosity (viscous diffusivity), p/p 

^ abscissa in auxiliary t-plane 

p fluid density 

$ dimensionless potential, <p/|v^|b 

(p potential function 

dimensionless stream function, ’/'/|v^|b 
i// stream function 

^ ^ r) 

V gradient, V = i — + j — 

9x 0y 

V dimensionless gradient, V = i — + j 

9X 0Y 

Subscripts: 

L along heated wall 



o reference value at 9 = tt/A 
s at free streamline 
w at wall 

°° condition of undisturbed incoming fluid 
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APPENDIX B 


MAPPING FOR JET FLOW IN A CORNER 
Mapping of Flow Between Physical and Potential Planes 

The flow being considered here is steady, incompressible, inviscid and irrotational. 
The region surrounding the jet is at constant pressure so that the outer jet boundary is a 
free streamline. The details of this type of flow, such as the velocity aloi^ the impinge- 
ment plate and the shape of the free streamlines, can be obtained by conformal mapping 
(ref. 4). A dimensionless system of variables is used as defined in appendix A. Addi- 
tional information on this type of solution is in reference 2. 

The flow is expressed in terms of a complex potential: 


W = $ + i’F 


(Bl) 


The derivative of in a direction gives the local velocity in that direction, and the 
derivatives of 'F are related to the derivatives of # by the Cauchy -Riemann relations. 

The derivative of the complex potential W provides the complex conjugate of the 
velocity in the flow region. 


= U - iV = C (B2) 

dZ 


By integrating, the relation between the physical plane and the potential plane is found as 


Z - / -idW + constant (B3) 

J ^ 

The known conditions of the flow can be used to determine the flow regions in the ^ and 
W coordinate planes. To integrate equation (B3) these ^ and W regions are each 
mapped conformally into an intermediate t plane. Once the W(t) and ^(t) conformal 
transformations are known, they can be siibstituted into equation (B3) and Z(t) obtained 
by integration. The integration constant is evaluated by matching a point in the t and 
Z planes. The conditions in the flow, such as the velocity, are then known as a function 
of position implicitly through the variable t. 
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Hodograph Plane and Transformation to Auxiliary t -Plane 


The flow region in the hodograph plane is shown in figure 15. The free streamline 
is in a region of constant pressime and hence has a constant velocity magnitude along it 
which is equal to unity by virtue of the dimensionless velocity at point 2. Hence, the 
free streamline transforms into part of a unit circle in the ^ -plane. The velocity is 
zero at point 4, and the velocity direction must be along the axis of symmetry 34 and 
along the plate 45. The hodograph is thus a sector of a unit circle. 

A convenient geometry for the intermediate t plane is a quarter circle as shown in 
figure 16. The mapping between the ?- and t-planes is given by (ref. 4) 

?(t) = (B4 


Potential Plane and Transformation to Auxiliary t- Plane 
From one of the Cauchy -Riemann equations and the definition of velocity potential 

^ - ^ = -Y (] 

ax 0Y 


Between points 3 and 2 in figure 3(b) V = -1 so that by integrating 


/•X(2) fX{2) 

4^(2) - ^^(S) = / -V dX = / 

^X(3) yx(3) 


dX = 1 


(B61 


The ^(3) is arbitrarily set equal to zero and then ^(2) = 1. The line 345 is thus the 
zero streamline, and the free streamline 12 is the unit streamline. The potential plane 
as shown in figure 4 is thus an infinitely long strip of imit width. 

The transformation of a unit strip (W -plane) into a quarter circle (t-plane) is given 
by 


W(t) 



(B7) 
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Integral to Determine Z(t) 


Noting that the condition Z = 0 when t = 0 (point 4 in figs. 3(b) and 16) yields 
zero integration constant, the relation in equation (B3) can be written as 


■/' 


1 dW 


'o ?(t) clt 


dt 


Differentiating the W(t) in equation (B7) gives 


^ a 

dt 


»(1 - t ^) 


Then 


Z = 


f 8t _ 8 „le(»/2) f dt 

Jg tW ,(1-A ” Jg 


Along the heated plate 45 the t is real and equal to ^ so that 


/ 


Z = 8 gi^(7T/2) 

" /o 




Z on 45 


The coordinate L aloi^ the plate is related to Z by 


Z = Le‘^^'/2) 


Hence the L along the plate is given by 


L = 


8 


/ 




’ /o /a - (^) 


L along 45 


To relate L to equation (B7) gives along 45 where t = 4 


(B8) 


(B9) 


(BIO) 


(BID 


(B12) 
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and since ^ is real this gives 


W = $ + i'Jf =f.ln| 

It 


k±£ 

,!-{■ 



7T 



This can be rearranged into the form 

$ = i tanh"^(4^) 0 :< ? s 1 

ir 


(B13) 


(B14) 


Equations (B12) and (B14) relate L to $ by use of the intermediate variable this 
relation will be needed to transform the wall temperatures obtained in the W-plane into 
the physical plane. 


Free Streamline 

To obtain the equations for the free streamline 12, it is noted that the streamline is 
along the unit circle in figure 16 so that 

t = (l)e^® 0 ^ 6 < - 

2 

Since the streamline starts and ends at infinity it is convenient to start the integration 
for Z in equation (B9) at a point (corresponding to e^) on the streamline, (The 

Xo, Yo will be determined later. ) Using equation (B9) and the fact that along the stream- 
line dt = ie^® de gives 


t ; - 8 Jl3iir/2) 
Z^__e 


/ 


le 


ie 


ei^6(i _ gi46) 


de 


where Z. is the value of Z at 6 = e 

° -i2fl 

numerator and denominator by e ^ 


o‘ 

to 


To simplify this expression multiply the 
obtain 
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I 







le 


-ie 


ei/30(e-i2e _ ^i2e 



e^^®(-2i sin 2e) 


-de 

) 


de 


# 



Separating real and imaginary parts 



-i|^e(i+0)-;8 Ij 

de 

sin 2e 


cos e(l + /3) - /3 - 

--L. Ll de 0 s e < 1 (Bisa) 

sin 2e 2 


Y-Yo 


4 

7T 


/ 

•'e=eQ=3T/4 



e(i + ;S) 


sin 2e 



0 e ^ - 
2 


(B15b) 


where the point has been chosen at e^ = tt/ 4 as a convenient point along the 

streamline. 

To determine X Y equation (B8) will be integrated from Z = 0, t = 0 to Z = Z 

u u ^ O’ 

t = e , It is convenient to carry out the integration along a line of constant 6 = 6. 

Then ° 


t = |t 
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and 


i0 

dt=d|t|e ° 




7T 



e'^°d|t| 


[t|=0 - |t|"e'"S 


Letting = tt/4 


Z = f ^1^1 

° • lu„|t|^(l..|t|4) 


Separating into real and imaginary parts 

Xq = ® cos (1 + )3)j F()3) 


Y„ =isin r?L(l + /3)1 f{)3) 
° 7T L4 J 


where 


F(/3) = 


f' d. 

»/ 7?^(1 + 7 ?^ 


(B16a) 

(B16b) 


Velocity Along Plate 

Along the plate 45 the velocity must be tangential to the surface so that 

^ = U - iV = |u[e 
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Also along 45, t = | so from equation (B4) 


Hence, along the plate 


^ ^ l+^g“i/3(ir/2) 

|U| =Ul = {*+« 


Then using equation (B 14) 


r 




U. 


tanh 


(l+/3)/2 


APPENDIX C 


MAPPING FOR INFINITE ROW OF JETS 

From the periodic nature of the flow in figure 2, it is necessary to consider only 
the portion of the flow between the two centerlines on either side of the y-axis. This 
portion is shown in the dimensionless physical plane in figure 5. 

The corresponding hodograph is shown in figure 17. The magnitude of the flow 
velocity is unity along the free streamline 345; hence this line transforms into a unit 
semicircle in the hodograph plane. Along the c^terlines 12 and 67 the U velocity 
component is zero, and along the heated plate 781 the V component is zero. The ve- 
locity is zero at the impingement and reverse stagnation points 7 and 1. 

In the potential plane (fig. 6) the free streamline becomes a horizontal line. Par- 
allel to it is the streamline 6712 formed by the centerlines and the heated plate. Thus 
the flow region transforms into a strip of unit width. The physical coordinates are re- 
lated to the hodograph and potential coordinates as previously derived in equation (B3). 

Z = J idW + constant (B3) 

To relate C and W in order to carry out the integration, the auxiliary u plane in 
figure 18 is used in which the flow region is mapped into a rectangle. The relation be- 
tween f and u can be obtained from results in reference 4 as 

g(u) = ^ ^ (Cl) 

k ' + dn u 


where 


k'^ = 1 - 


The mapping between the W and u planes is derived by use of the intermediate T 
plane in figure 19. The W is related to T by 


W = i In 





(C2) 


26 



The transformation between T and u is by means of an elliptic function; 


I 

3 


T = sn u 

Equation (B3) is written in terms of u as 


Z = / J- 

J ?(u) du 


du + const 


Substitute sn u for T in equation (C2) and differentiate to obtain 


dW 

du 


2 cn u dn u 
± - sn^u 


(C3) 


Then 


Z 


2 f k* + dnu kcnudnu ^^ ^ 
nj kcnu i_k2gA 


Using the identity 


1 - 


,22 .2 

k sn u = dn u. 


reduces this to 


Z 


2 

n 



+ const 


Integrating yields 


Z 



(C4) 


where the integration constant has been set equal to zero by use of the condition at 
point 8 that Z = 0 when u = 0. The spacing between the jets is related to K (and hence 
to k) by the fact that at u = K, Z = S/2b. Then from eqviation (C4), 


2b 



or 



(C5) 
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With K determined from a specified S/2b, the k is also known from the definition of 
the complete elliptic integral, 


K = 



d^ 

Vl - k^ sin^ <p 


0 < k^ ^ 1 


and the k* is found from 


k’ = Vl -k^ 

Since the minimum value of the elliptic function K is ir/2, equation (C5) shows that the 
minimum spacing consistent with potential flow is S/2b = 2. As the spacii^ is de- 
creased toward this value, it will be shown that the free streamline moves out toward 
infinite distance above the heated plate; that is, the space between the jets completely 
fills with fluid. 


Coordinates of Free Streamline 

From figure 18 the free streamline corresponds to u = ^ + iK*. Substituting into 
equation (C4) for Z gives 


Z(^ + iK') - i 


■'(— y 

Vdn kJ 


tan"^f-2i-!-]+ | + iK' 


which is equivalent to 


X + iY|^=£ 
7T 


llnf^ » . l -t- ^ -V|H-iK' 
2 Vdn f -k’/ 


The X,Y coordinates along the free streamline are then 


X = £^ -K<4<K 

IT 


Y = ^ [k’ + llnf-— ^ ^ \ 

7T L 2 Vdn 4 - k’/_ 


-K s ^ < K 


(C6a) 


(C6b) 
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Tangential Velocity and Relation Between # and X Along Heated Wall 


The solution of the energy equation involves the tangential velocity along the heated 
wall. Also the # along the heated plate in the potential plane has to be related to the 
X location in the physical plane. 

Along the heated plate -S/2h <X s S/2b, Y = 0, the u = ^ with the range 
-K < I < K. Then along the plate, equation (C4) becomes 

X = i (C7) 

Along the heated plate = 0 so that W = $. Also u = | so that from eqtiation (C3) 
T = sn Then usir^ equation (C2) 

4. (C8) 

V \1 - k sn 

Equations (C7) and (C8) relate the # and X throt^h the use of the intermediate varia- 
ble At point 1, ^ = K so that sn ^ = 1 and then 

^ (point 1) = 1 In / (C9a) 

V Vl-k/ 

At point 7, I = -K so that sn ^ = -1 and then 

$ (point 7) = 1 In / ^ (point 1 (C9b) 

IT Vl + k/ 

From equation (Cl), since the V component is zero along the plate, the U com- 
ponent along the plate becomes 

U=-A£!LL_ (CIO) 

k’ + dn ? 

Equation (C 8) is solved for sn | to yield 


sn 4 = — 


1 1 - e“ 


k 


= -itanh 
k 
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Then using the identities cn ^ = Vl - sn^^ and dn ^ = Vl - sn^^, the expression 

for U along the plate becomes 



(Cll) 
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APPENDIX D 


SOLUTION FOR WALL TEMPERATURE AT LARGE I TOR CORNER FLOW 


Let the coordinate system be along and normal to the impingement plate as 

shown in figure 20. At large I the flow distribution has become uniform with velocity 
v^. The convective heat-transfer equation in the region of large S. is 


VoopC 


P0£ 



(Dl) 


hi this region the temperature distribution has become fully developed, that is, the shape 
of the temperature profile across the n^ direction is no longer changing with increasing 
£. As a result of the uniform heat addition at the wall, however, the temperatures 
across the width of the flow are all increasing linearly with f , The temperature distri- 
bution for large £ must then have the form 


t(f,n^ = + f(n^) 


(D2) 


This form of the solution is substituted into the energy equation (Dl) to yield 


d2f 


VoopCp^’i - « — 


(D3) 


To obtain f integrate once to give 


df 

^«w 


VoopC 


- Vw + ^2 


To satisfy the condition that d£/dn = 0 at the Insulated boundary n„^ = b, the 
V pC '^’i 

= - °° ^ _P- l b. 

K 


A second integration then yields 




(D4) 
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To obtain and ‘ifg an overall heat balance is used from the incoming jet at 
large y in figure 20 where the fluid is all at t^, to a location at large £ where the 
mixed mean flmd temperature is t(£). The heat balance yields 

= bv^pCp[t(£) - t J (D5) 

In the fully developed region from equations (D5) and (D2), 


= __2L_= 

9£ 0£ “bv^pCp ^ 


(D6) 


so that the f in equation (D4) becomes 




The temperature distribution in equation (D2) becomes 


t(£, V = 


q f 

^w 




Kb \ 2 




+ 


(D7) 


For a uniform velocity distribution the mixed mean fluid temperature t(£ ) is given by 




(D8) 


Then, substituting equation (D7) into (D8) gives 



t(£) = 


V 

toVeopCp 


3k 




(D9) 
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Equating this to t(£) from equation (D5) gives 


bVeopCp^ bv^pCp 


3/c 


+ <5f3 


Hence, 


«'3=*oo + 


3k 


(DIO) 


The fluid temperature distribution in equation (D7) then becomes 


t(^, n^) = 




bv^pC Kb V 2 


-fhr-»"w 


+ + 


q b 
3fc 


To obtain the wall temperature, this is evalxiated at = 0 

V<«“: 




bv^C 

hi dimensionless form the final result is 


3k 


(Dll) 


(D12) 


T„(l.) - T„ = A L a (DIS) 

which is valid for large L. Figure 8 indicates that, depending on the plate angle, L 
can be in the range of about 2 to 10 to obtain reasonable agreement with the complete 
solution. 
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Heated wall 


(a) Flow into internal corner (y positive). (b) Flow over external corner (y negative) 

Figure 1. - I nviscid-ir rotational slot jet impinging in corner region formed at intersection of two 
uniformly heated plates. 



Figure 2. - Impingement of infinite row of liquid-metal slot jets against 
uniformly heated surface. 






(a) Physical plane. (b) Dimensionless physical plane. 

Figure 3. - Jet region in physical plane for impingement in corner. 



Figure 4. - Jet region and boundary conditions for flow in corner region mapped into potential plane. 
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Figure 5. - Jet region in dimensionless physical 
plane for infinite row of impinging jets. 




5 


r\p = i 


( 0 , 1 } 


3 



Figure 6. - Jet region and boundary conditions for infinite row of jets mapped into potential plane. 
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Dimensionless coordinate, X = x/b 

(a) Impingement into internal corner {positive y). (b) Impingement over an external corner (r\egative y). 

Figure 7. - Free streamlines for jet impinging In corner region. 
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Dimensionless wall temperature distribution 





Dimensionless wall temperature distribution, 


11 



(b) Jet Peclet number, 10. 
Figures. - Continued. 
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Dimensionless wall temperature distribution 





Dimensionless wall temperature distribution 
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(a) Jet Peclet number, 5. (b) Jet Peclet number, 10. 

gure 9. - Effect of angle of impingement plate on dimensionless wall temperature distribution for fixed value of iet Peclet number at plate angles 
from -90 to 0° 




Dimensionless distance along heated plate, L = l/b 


(c)Jet Peclet number, 20. (d) Jet Peclet number, 50. 


Figure 9. - Concluded. 




stagnation heat-transfer coefficient, j- ; _ Dimensioniess wail temperature at stagnation 


20 , 
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Figure 10. - Dimensionless wall temperatures at stagnation (apex of flat plates) as function of plate angle for 
various jet Peclet numbers. 
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Figure 11. - Stagnation heat-transfer coefficient for plate at angle, relative to coefficient at stagnation point of 
flat plate. 
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Figure 19. - Auxiliary T-plane. 
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Figure 20. - Geometry for derivation of heat-transfer behavior at large i. 
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